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@ Proc6d6 de correction spatial e d'Image num6rique. 

(57) Un proc^dd de correction spatlale d* image 
num^rique compos^e de rang^es de pixels qui 
sont des lignes ou des colonnes, comprend la 
comparaison (E20) des rang^ deux k deux 
pour determiner des premiers d^caiages entre 
rang^es exprim6s en nombre entier de pixel, la 
determination de seconds d^calages entre ran- 
ges exprim6s en nombre non entier de pixel 
puis la determination (E40) pour chaque rangee 
d'un troisieme d^calage par rapport k une ran- 
g6e de reference, et enfin la correction de 
chaque rangee par un decalage des pixels de 
cette rangee par rapport d la reference egal ^ la 
partie entiere du troisieme decalage de cette 
rangee et par un lissage spatial de ses pixels 
avec des coeffidents fonction du troisieme 
decalage. Ce precede ameiiore la quallte de 
rimage lorsque celle-ci a ete prise par un aero- 
nef soum'» d du roulis et/ou de la derive. 
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La pr6sente invention concerne d'une mani^re 
g6n§rale le traitement d'innage num6rique pour amd- 
liorer I'aspect visuel de d'lmages et ainsi en fadliter 
Texploftation ult6rieure. 

L'inventlon a trait plus partieuliferement d un pro- s 
c6d6 de correction spatiale d'iniages num6riques 
capt6es par des disppsitlfs enibarqufe dans des en- 
gins mobiles tels qu'avions ou navires. Un tel dispo- 
sitif "photographle' une mage ligne par ligne et re- 
constitue riniage par la juxtaposition des lignes. Or io 
rengin mobile subit au cours de son d^placement des 
mouvements perturbants tels que roulis ou derive (la- 
cet) qui provoquent des distorsions de IMmage. 

Si le dispositif embarqu6 est fixe par rapport d 
rengin nnobile, un gyrom6tre de mesure de vitesse is 
angulaire peut 6tre associd au dispositif de manite 
d n^esurer en temps r§el Tangle des mouvements per- 
turbants par rapport ^ la trajectoire de I'engin et d 
corriger en consequence la distorsion de Tlmage. 

Alternativement, le dispositif embarqu6 peut Stre 20 
place sur une plate-forme stabilisee qui corrige tout 
mouvement perturbant Le dispositif embarqu6 ne 
subit alors plus de mouvement perturbant 

Ces materiels de connection, tels que gyromdtre 
de mesure de vitesse angulaire ou plate-forme stabi- 25 
Iis6e, sont tr6s coOteux et ne corrigent pas des dis- 
torsions suppiementaires, sous fomne d'instabilit6 
dans le sens ligne de Timage. provoqu6es par du bruit 
sur un signal de synchronisation ("gigue") lors de la 
transmission des images par exemple de I'engin mo- 30 
bile, tel qu'avion ou navire, vers une station terrestre 
de traitement et diffusion ou reproduction d'image. 

La pr6sente invention vise d fournir un proc6d6 
de traitement spatial d'image num6rique qui am6llore 
raspect visuel d'une image de manl6re 6conomique, 35 
en corrigeant non seulement des distorsions d'image 
dues ^ des mouvements perturbants, mais dgale- 
ment des distorsions dues d des d^fauts de synchro- 
nisation. 

A cette fin, un prooldl del^orreSion sp^^ 40 
d'image num6rique compos6e d'un ensemble de ran- 
g6es de pixels qui sont des lignes ou des colonnes de 
Timage, est caracterise en ce qu'il comprend 

la comparaison des rang6es deux ddeu^ 
determiner des premiers d6caJages entre rang^es 4S 
exprimes en nombre entler de pixel, 

la determination de seconds decalages entre 
rangees exprimes en nombre non entler de pbcel. en 
fonctlon des premiers decalages, 

la determination pour chaque rangee d'un troi- so 
sieme decalage de ladite rangee par rapport d une 
rangee de reference en fonctlon des seconds deca- 
lages, et 

la correction de chaque rangee par un decala- 
ge des pbcels de ladite chaque rangee par rapport d 55 
la rangee de reference egal d la partie entiere du trol- 
sieme decalage pour ladite chaque rangee, et par un 
lissage spatial de pixels de ladite chaque rangee avec 



des coefficients de lissage fonctlon du troisieme de- 
calage pour iadlte chaque rangee. 

La comparaison des rangees deux e deux 
comprend pour deux rangees donnees, le calcul de 
sommes de correlation entre les deux rangees deca- 
I6es respectivement lorsque I'une est soumise par 
rapport d I'autre e des decalages predetermines ex- 
primes en nombre entier de pbcel, la recherche d'une 
somme de correlation minimale parmi les sommes de 
correlation calcuiees afin d'en deduire ledit premier 
decalage entre les deux rangees exprime en nombre 
entier de pixel. 

Une somme de correlation entre deux rangees 
decaiees d'un decalage exprime en nombre entler de 
pbcel est egale e la somme des valeurs absolues des 
differences des valeurs des pixels d'une rangee et 
des valeurs des pbcels decaies dudit decalage de I'au- 
tre rangee. 

De preference, la detemnination du second deca- 
lage entre deux rangees exprimee en nombre non en- 
tier de pixel est af f inee par une interpolation des som- 
mes de correlation fonctlon des decalages predeter- 
mines entre les deux rangees, par une fonctlon pdy- 
ndmiale d'ordre superleur ou egal e 2 autour du pre- 
mier decalage, le second decalage etant I'absclsse 
d'un minimum de la fonctlon d'interpdation. Par 
exemple, la fonctlon polyndmiale peut-etre une fonc- 
tlon parabollque dont le sommet oonstltue ledit mini- 
mum. 

La detenminatlon pour une rangee donnee d'un 
troisieme decalage de ladite rangee par rapport d une 
rangee de reference est r6allsee par la somme des 
seconds decalages des rang6es comprises entre la 
rangee de reference et la rangee donnee. 

Selon une autre caracteristlque, les coefficients 
de lissage spatial d'une rangee sont tels que leur 
somme est egale e 1. Par exemple, les coefficients 
de lissage spatial d'une rang6e sont tels qu'un pre- 
mier coefficient est egal ^1 la partie dedmale du troi- 
sieme decalage de la rangee et qu'un second coeffi- 
cient est egal ou compiementaire d 1 du premier coef- 
ficient 

De preference, les seconds decalages entre ran- 
gees sont remplaces par des decalages coniges 
fonctlons des seconds decalages et tels que leur va- 
leur moyenne est nulle. Cette etape est fbndee sur la 
distinction entre les decalages e conriger dus d des 
perturbations telles que roulis et des "decalages na- 
turels" dans I'lmage, tels qu'une route ou un cours 
d'eau non paralieie e la trajectoire de I'engin. 

En variante, lorsque les decalages e corriger sont 
dus d de la derive, les seconds decalages entre ran- 
gees sont remplaces par des decalages corrig6s res- 
pectivement 6gaux d la moyenne des seconds deca- 
lages. 

D'autres caracteristiques et avantages de la pre- 
sente Invention apparallront plus dairement k la lec- 
ture de la description survante de plusieure realisa- 
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tions pr^f^r^es de rinvention. en r6f6rence aux des- 
sins annexes con-espondants dans lesquels : 

- la figure 1 est une vue sch^matique en pers- 
pective d'un a^ronef 6qu(p6 d'un disposltif de 
prise de vue; 

- la figure 2 est une representation sch^matique 
d'un premier type de deformation d'Image; 

• la figure 3 est une representation schematique 
d'un secx)nd type de deformation d'Image; 

- la figure 4 est un bloc-diagramme des etapes 
principales d'un precede de correction d*image 
selon une premiere realisation; 

- la figure 5 est un atgorithme de determination 
de premier et second decalages entre lignes 
selon une premiere realisation; 

- la figure 6 est un algorlthme de correction de 
second decalage entre lignes selon la premie- 
re realisation; 

- la figure 7 est un algorlthme de determination 
de rimage corrigee selon la premiere realisa- 
tion; 

- la figure 8 est un bloc-diagramme des etapes 
principales d'un precede de correction d'Image 
selon une secende realisation; 

- la figure 9 est un algorlthme de determination 
de premier et second decalages entre lignes 
selon la secende realisation; 

- la figure 10 est un algorlthme de correction de 
second decalage entre lignes selon la secende 
realisation; et 

- la figure 11 est un algorlthme de determination 
de rimage corrigee selon la secende realisa- 
tion. 

En reference k la figure 1 , dans un aerenef 1, tel 
qu'un avion ou un drone, est fbce un meyen de prise 
de vue 2 comprenant essentiellement un nnoyen de 
detection de preference sensible aux inf rareuges, un 
analyseur mono-ligne et un moyen de memoire et/ou 
un meyen de transmission de dennees. 

L'aeronef 1 survele une bande rectiligne de 
terrain 3 de maniere e phetographlercette bande. La 
prise de vue est effectuee par un balayage eiectronl- 
que llgne par ligne. Une llgne L est une "bande eie- 
mentaire" de la bande 3, transversals e la bande 3. 
et la Juxtaposition d'un ensemble de lignes successl- 
ves forme une image d'une portion de bande 3, le de- 
placement de I'aerenef le long de la bande fournrs- 
sant le balayage d'Image (frame). 

Une ligne captee par le moyen de detection subit 
notamment une conversion anaioglque/numerlque en 
entree de I'analyseur. Une ligne est un ensemble de 
M pixels et un groupe de N lignes est une matrice de 
N X M pixels, foment une image num6rique e N lignes 
et M colonnes. Une ligne ou une celonne de Timage 
constitue une rangee de pixels dans la matrice, et par 
extension une rangee de I'image au sens de la pre- 
sente invention. Une ligne numerique est memorisee 
par le moyen de memoire. En variante, les niveaux de 



lumlnosite des pixels de la ligne sent emis par le 
moyen de transmission vers un recepteur associe e 
un moyen de memorisation et de traitement de don- 
nees. Le terme "lumlnosite" se rapporte e la tuminan- 

5 ce pour une Image en neir et blanc, ou bien e Tune des 
cempesantes prinoaires rouge, verte et bleue ou e Tun 
des signaux de luminance et de chrominance pour 
une image en couleur. 

Le dispositif de prise de vue et en particulier le 

10 moyen de detection optique sent fixes par rapport e 
Taeronef. Pour obtenir des Images reproduisant le 
plus f ideiement possible la scene effective de la ban- 
de 3, I'aerenef vole de preference en ligne dreite et k 
altitude constants. Cependant, en raisen des oondl- 

15 tlons atmespheriques, i'aeronef subit en cours de vol 
des mouvements perturbants notamment de roulls, 
c'est-e-dire des petites rotations alternees R autour 
de i'axe longitudinal de I'aerenef, et de derive, c'est- 
e-dlre une rotation D autour d'un axe vertical, ou de 

20 lacet, c'est-e-dire des petites rotations altern6es au- 
tour d'un axe vertical. 

Dans la figure 2 est mentree schematlquement 
une image perturbee uniquement parle roulls de I'ae- 
renef. La forme geemetrique reelle de la scene sur la 

25 . portion de bande de terrain photographiee est un rec- 
tangle fence, sur un fend dair. L'image fburnie par 
I'analyseur e bord de I'aerenef soumis au reulis pre- 
sente une forme geemetrique defermee. 

L'image est composes de N lignes Li e L^ et de 

30 M oolennes Ci e Cm. dont une llgne et une celonne 
Cm sent representees e titre d'exemple, avec 1 ^ n ^ 
N et 1^ m ^ M. Le reulis proveque une deformation 
de I'image dans le sens llgne, c'est-e-dire les lignes 
sent decaiees les unes par rapport aux autres. 

35 En pratique, une Image comprend par exemple 
1024 lignes et 1024 colonnes. 

En reference e la figure 3, une image schemati- 
que est perturbee uniquement par des derives de 
I'aerenef principalement suivant un meme sens par 

40 rapport e la trajectoire de I'aerenef. La forme geeme- 
trique reelle est e neuveau un rectangle fence sur 
fend dair. L'image feurnie par I'analyseur e bord de 
I'aerenef soumis k la derive presents une forme gee- 
metrique defermee resultant de la transformation du 

45 rectangle d'erigine senslblenoent en lesange. 

Une ligne Ln parmi N lignes et une celonne C^ 
parmi M colonnes sent egalement representees dans 
la figure 3. La derive proveque une deformation dans 
le sens celonne, c'est-e-dire les colonnes sent deca- 

50 lees les unes par rapport aux autres. 

En pratique, les deformations dues aux reulis et 
les deformations dues k la derive se combinent aiea- 
toirement. Dans le cas d'images reelles, porteuses en 
general d'un grand nombre de details, ces deux types 

55 de defemiation rendent diff idles la lecture et I'inter- 
pretatlen des images parun operateur. De plus, d'au- 
tres deformations alterant une image. En particulier, 
dans le cas de la transmission d'images entre I'aero- 
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nef et un moyen de reception, la synchronisation des 
donn6es est r6alis6e par une r6f6rence de seuil sur 
unfrontmontantd'un s^nal de synchronisation. Lors- 
que le front nriontant est bruit6. la position du seuil sur 
le front montant subit des variations provoquant un 
ph6nom6ne de gigue Oitter). c'est-S-dire d'instabilit* 
dans le sens des llgnes de rimage. Cette instability 
provoque des d^calages entre lignes analogues aux 
dWormations dues au roulis. mais de plus foibles am- 
plitude. 

En r6f6rence d la figure 4. un proc«d6 de 
correction d'bnage destln6 plus particuliferement d 
une Image d6fbrm6e par du roulis est iliustr^ sous 
la forme d'un algorithme compos* de anq 6tapes 
pnncipales E10 k E50 et programme dans une nZ 
moire de calculateur. La correction d'image est rtell- 
see dans ra^ronef, ou en variante au sol. Dans hxis 
les cas. comme ii a d«jd 6t6 pr6cis6, une image est 
une matrice de N x M pbcels m6moris6e. Le proc6d6 
comprend une premifere 6tape E10 de lecture d'une 

•magequiconsisteilirelesniveauxdeluminosMdes 
N X M pbcels de I'imaga i N lignes et M colonnes dans 
une mdmoire d'image. 

L'«ape E20 determine et m6morise les d6cala- 
gw des lignes prises deux d deux successivement 
L6tape E20 est compos6e d'un ensemble d'^tapes 
E21 a E27 qui seront dtoites en detail en r«f«ience 
a la figure 5. 

L'«apeE30traiteles d6calages entre lignes pour 
S6parer les d6calages dus au roulis, de tout d6calage 
non h6 d la prise de vue. mais existant r6ellement sur 
la portion de terrain photographic, par exemple une 
route formant un angle avec les boids de I'image. ou 
un cours d'eau qui serpente. L'tepe E30 est compo- 
st d un ensemble d'dtapes E31 d E34quiserontd*. 
cntes en detail en reference d la figure 6. 

L'6tape E40 determine I'image corrig6e k partir 
des d6calages entre lignes d6termin6s d I'dtape E30 
En variante, l'6tape E40 determine I'image corrigte |) 
partir des d6calages entre llgniis d6termln§s d i'6tape 
EZO; cette variante ne comprend pas l'6tape E30 

Dans lous les cas. I'6tape E40 se compose d'un 
ensemble d'6tapes E41 d E46 qui seront d6crites en 
detail en r6f6rence k la figure 7. 

La nouvelle image corrig^e des ef fets du roulis 
est memorise h retape E50, et est. de pr6f6rence. 
6galement aff ich6e k l'6cran d'une console de visua- 
lisatlon associde au calculateur. 

En variante, I'image corrig6e k l'6tape E40 peut 
etre trait^e k nouveau k partir de I'^tape E10 Hu- 
sieurs cydes de correction successifs am^liorent 
plus I image qu'un cyde unique. Le nombre de cycles 
de correction est soil fixe, d6termln6 une fois pour 
toute. par exemple 6gal in. 2 ou 3. soil choisi par 
lop^rateur, soil encore d6termin6 dynamiquement 
par I algorithme sur un crit^e de quality de I'image 
comgee. ' 

Comme illustrte d la figure 5. la determination 



dM dtealages entre lignes comprend un ensemble 
d 6tapes E21 k E27. Une ligne quelconque L„ de 
I image est compares k la ligne adjacente , pr6c6. 
demment m6moris6e. L'6tape E21 initialise la deter- 

5 mination k la ligne Lj qui est compai^e k la ligne L, 
Dans la suite, la ligne L, n'est pas affectee par to 
correction d'image et est une ligne de reference. En 
vanante. Fordre des calculs est imrers6, la ligne de re- 

« . ^ ®' «»"Paraison commence 

10 alors k la ligne U qui est compar6e k la ligne U 
Selon une autre variante. la ligne de reference est 
une ligne quelconque de I'image. 

L'etape E22 Initiailse un decalage en nombre en- 
tier de pixel Dec k une valeur-Dmin. Le decalage Dec 

IS est le decalage de comparaison de la ligne L, avec la 
ligne L,^,. 

Pour une valeurde decalage Dec donnee. l'etape 
E23 compare la ligne I, k la ligne „ en calculant 
une somme de comeiatfon So„ entre les deux lignes 
20 pour le decalage Dec. La somme Soee est egale a: 

m^K-Bord 



oO: 



30 



35 



40 



- P'Xd- i,in est le niveau de luminosite du pbcel au 
croisementde la ligne L^ , etde la colonne C„. 

i^l^nttoae est le niveau de luminosite du pixel 
au croisement de la ligne L, et de ia colonne 

- Bord est un entier, strictement Inferieur k M 
6gal k un nombre de pixels situes sur chacun 
des deux bonis de i'image. et non consideres 
dans la comparaison; k titre d'exemple, pour 
des lignes de 1024 pixels, les boiris non consi- 
deres sont d'une largeur comprise entre 50 et 
100 pixels. 

La valeur de est memorises, puis le decala- 
ge de comparaison Dec est incremente de 1 

AI'etapeE24.le decalage Dec est compare e une 
valeur predeterminee Dmax qui est la valeur maxima- 

ledudecalageDec. Si la valeurdu decalage Dec est 
4S infeneure ou egale k la valeur Dmax, l'etape E23 est 
k nouveau derouiee. Sinon. toutes les sommes de 
correlation So., pour Dec variant de -Dmin k Dmax 
sont detemiinees, c'est-e-dire des sommes au nom- 
bre de (DmIn + Dmax + 1). et I'algorithme passe k 
so I etape E25. 

Les bornes de variation -Dmin et Dmax du deca- 
lage Dec sont k titre d'exemple -25 et +25 pbcels 

L'etape E25 recherche la somme minimale Sn.^ 
panni les sommes determinees precedemment. ainsi 
« que le decalage Decmin associe. Pour trouver la 
somme minimale Sd«^ parmi (Dmin + Dmax + 1) 
sommes, les valeurs adjacentes kS^^sont compa- 
rees deux k deux. La somme minimale So^^ est la 
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premiere somme trouv6e dont les deux sommes ad- 
jacentes pr6c6dente et suivante lui sont sup^rieures 
lorsque Dec crolt de - Dmin d Dmax. En variante, lou- 
tes les sommes Soec sont anatys6es, en comparant 
les sommes deux d deux, et en gardant comme r6sul- 
tat de chaque comparaison la somme Sdm Plus fai- 
ble. Selon une autre variante, la somme minimale 
Soecmin 6st la premiere somme trouv^e dont les deux 
sommes adjacentes pr^c^dente et suivante lul sont 
sup^rieures de valeurs de difference sup6rieures^ 
un seuil predetermine. 

Dans tous les cas, d la soimme minimale trouvee 
^Decmin associe le decalage de la ligne L„ par rap- 
port d la ligne U>. i exprlme en nombre entler de pixel 
Decmin. Or, un decalage exprime en nombre entier 
de pixel ne sert d corrlger une Image que d'un deca- 
lage relatlvement eieve d'une ligne d Tautre. En pra- 
tique, le decalage d'une ligne d i'autre survenu au 
cours de la prise de vue est f requemment inferieur au 
pbcei; rinventlon prevoitalorsd'af finer le calcui du de- 
calage. 

L'etape E26 determine le decalage Dn de ta ligne 
Ln par rapport ^ ia ligne Ln. i, avec une precision infe- 
rieure au pixel; d'une maniere generate, Dn est un 
nombre non entier de pbcel. Le decalage Dn est deter- 
mine en fonctlon des decalages Decmin, (Decmin-1) et 
(Decmin 1) ainsl que des sommes de correlation as- 
sociees d chacun de ces decalages Doecmin , Soeonirv i 
^oecmift^i- Dans un repere ayant les decalages pour 
abscisses x et les sommes de conreiation pour ordon- 
nees y, ces trols couples de decalage et somme de- 
finissent trois points. En recherchant la parabole 
d'equation : 

y = ax2 + bx + c 
passant par ces trois points, est deduit le som- 
met de cette parabole qui a pour ordonnee la somme 
minimum affinee et pour abscissa le decalage plus 
precis Dn de la ligne par rapport d la ligne l-n. i. Le 
decalage D„ resulte de la resolution du systeme ^ 
trois equations de parabole correspondent aux trois 
points af in d'en deduire les coefficients a, b et a Ce 
decalage est ainsi egal ^ : 

- b/(2a) = Decmin + 0,5x(SDocmin + 1 - SDocmm- i) 

/ (2 X Soecmln * Soecmln - 1 ' Soecmln + l) 

La valeur de Dn est ensuite menfH)risee et le pa- 
rametre n est incremente de 1. 

En variante, le decalage Dn peut etre determine 
comme Tabscisse du minimum de la parabole pas- 
sant par les trois points (Decmin, Soecmin)* (Decmin -i- 
k. Soecmin ♦ k). (Docmin-k, Soecmin- 0 0^1 k est un entier, 
valant 2 ou 3 par exemple. Selon une autre variante, 
la parabole est remplacee par une fonction polynd- 
miale d'ordre plus eieve que deux. 

L'etape E27 compare la valeur de n avec le nom- 
bre total de lignes N. Si n est inferieur ou egal e N, il 
reste au moins une ligne pour laquelle le decalage Dn 
n'a pas ete detennine et le calcui reprend pour une 



nouveile ligne d partir de retape E22. Sinon, tous les 
decalages D2 ^ D^ sont determines et I'algorithme 
passe ^ retape E31 qui est decrite dans la suite. 
En variante, une ligne quelconque L.n n'est pas 

5 comparee S la ligne adjacente t, mais d la ligne L^. 1, 
oil i est un entier superieur e 1 , puis le parametre n est 
incremente de I. Atitre d'exemple, rentier I vaut 2 et 
seules les lignes de rang pair sont traitees. Les lignes 
de rang impair sont decaiees (etape E40) d'un deca- 

10 lage egal au decalage calcuie pour la ligne de rang 
pair precedente. Le nombre de calculs est reduit et ia 
Vitesse d'execution est augmentee. 

En reference d la figure 6, le traltement des de- 
calages entre lignes comprend quatre etapes prin- 
ts cipales E31 e E34. Ce traitement a pour but de corri- 
ger les valeurs des decalages Dn, pour n variant de 2 
d N, determinees k retape E26. En effet, les decala- 
ges Dn dependent des decalages dus au roulis, mais 
egalement de "decalages" representee par des traits 

20 inclines ou brises et existent reellement sur la bande 
de tenrain consideree, par exemple une route non pa- 
ralieie aux bords de la bande de tenrain. La correction 
des decalages repose sur le constat que les varia- 
tions de I'orientation de I'aeronef par rapport d sa tra- 

25 jectob-e dues au roulis sont nulles statistiquement 
Par consequent, les decalages dus au roulis sont nuls 
en moyenne sur un grand nombre de lignes, 1 024 par 
exemple, tandis que les 'decalages existent reelle- 
ment" sur le terrain ne sont pas a priori nuls en moyen- 

30 ne. La correction a pour r6sultat des valeurs de de- 
calage corrigees, dont la moyenne est nulle. La 
correction est du premier ordre sur les decalages; en 
variante, la correction est d*ordre plus eieve. 

La correction des decalages est telle qu'un en- 

35 semble predetermine de lignes de I'image n'est pas 
modif le lors de la determination de I'image corrigee, 
qui sera decrite en detail dans la suite. Atitre d'exemv 
pie, si rentier N est pair, cet ensemble predetermine 
de lignes comprend les lignes Li, Ln^ et l^. En varian- 

40 te, I'ensemble comprend un nombre plus important 
de lignes, par exemple 4,5 ou 6, qui sont reguliere- 
ment reparties ou non dans I'image. Selon une autre 
variante, le nombre et/ou la position dans I'image des 
lignes de cet ensemble sont chotsis par I'operateur. 

45 L'etape E31 determine trols coefficients St, S2 et 
P utilises pour la correction des decalages. Le coef- 
ficient Si est la somme de tous les decalages Dn, pour 
2^n^N. Le coefficient S2 est la somme des decala- 
ges Dn, pour 2^n^N/2, c'est-d-dire pour la demi- 

50 image entre les lignes Li et Lf^. Le coefficient P re- 
presente un coefficient egal d 4(Si-2S2)/N2 relatif e la 
variation du decalage des lignes par rapport d la pre- 
miere ligne de reference entre les lignes llmites Li et 
Lnq. et Lkq et Ln- 

55 L'etape E32 initialise la determination de la 
correction de decalage e la deuxieme ligne puisque, 
comme il a dej^ ete precise, la premiere ligne de 
I'image n'est pas modif lee lors de la determination de 

5 
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rimage corrig6e. L'6tape E33 calcule un coeff ident 
correcteur Cor^ qui est 6gal d S1/N + P x (n - N/2) af in 
de ramener la variation du d§calage des llgnes par 
rapport k la premfere llgne de r6f6renGe le long d'une 
droite, Le coef f ident Cor„ est ensuite soustrait du d§- s 
calage D„. ce qui prodult un d^calage corr^d Dc„. Le 
paramdtre n est mcr6nr)ent6 de 1. 

Wtape E34 compare le paramMre n d la valeur 
N; si n est inf6rieur ou 6gal d N, l»6tape E33 est d nou- 
veau d^oul6e pour cette nouvelle valeur de n. Sinon w 
tous les d6calages corrigds Dc„. pour 2^n^N. son! 
d6tennin6s et ralgorithme passe d T^tape E41 qui est 
d6taill^ dans la suite. 

En r6f6rence d la figure 7. la determination de 
rimage corrigde comprend un ensemt)l8 d'6tapes 15 
E41 d E46 et est effectu6e llgne par ligne. 

L'6tape E41 initialise la correction de llmage d la 
deuxifeme ligne L2 de rimage. Ar6tape E42 est d6ter- 
min6 un d6calage Decain de la ligne U qui est 6gal k 
la somme des d6calages Dc„ pour ^l^n, et qui d6. 20 
f init le d6calage de la ligne U par rapport d la ligne 
de r6f6rence L,. En variante. Decal„ est la somme des 
d6calages D, pour 2sslsn. calcuI6s d r6tape E26 
Dans tous les cas, le ddcalage Decain est en g^n^ral 
un nombre non entier, dont la partie entidre est not^e 25 
En. La partie d6dmale de Decal„ est not^e k„, et le 
compl6mentaire d 1 de k„i est notd k„2. Wtape E43 
initialise d 1 un param^tre m. repr6sentatlf du rang 
d'un pixel dans la llgne L„. 

L'6tape E44 determine le niveau de luminance 30 
conng6 Pc^ du pixel de rang m dans la ligne L„ par 
un d6calage de pixel 6gal d la partie enti^re de Decal„ 
puis par un lissage spatial de valeurs de pixel de la 
ligne. Le lissage rdsulte de la contribution de deux 
pbtels enregistrds adjacents de rangs m + En et m + 35 
En + 1 dans le pixel conig6. Le niveau corrig6 Pc,^„ 
est ainsi la somme des niveaux de pixel PIX„,,„+En et 
P'^m*En*i respectivement pond6r6s par les coeffi- 
cients kni et Kz, Le niveau corrlg6 Pc„,„ est m6moris6 
et le paramdtre m est Ihcr6ment6 de 1. En variante 40 
les coefficients k^^ et k„2 sont determines different 
ment. et sont par example 6gaux ^i 0 et 1 respective- 
ment, ou ^ 0,5 et 0,5. Selon une autre variante, le ni- 
veau corrige PCn^ est une somme ponder6e de plus 
de deux niveaux de pbcel. par example PIX„^^,. 45 
'^l^mEn et PIX^En*,. Dans tous les cas, les coeffi- 
cients de ponderation sont de preference tels que leur 
somme est egale dl.de sorts qu'il n'y art pas de mo- 
dification de la valeur moyenne des pbcels de la ligne. 

L'etape E45 verif ie si le parametre m est inf6rieur 50 
ou egal au nombre M de cdonne de I'image et done 
au nombre de pixels de la ligne L„. Si m est inferieur 
ou egal d M, il reste au moins un pixel dont la valeur 
conigee n'a pas ete determines, et retape E44 est d 
nouveau parcourue. Lorsque tous les pbcels de la ii- 55 
gne L„ sont corrig6s, alors le parametre n est incre- 
mente de 1 pour passer d la ligne suivante. Letape 
E46 v6rif ie si le parametre n est inferieur ou egal e N. 



En cas de reponse positive, I'algorithme deroule les 
etapes E42 e E45 pour une nouvelle ligne. Lorsque 
toutes les lignes sont comgees, ralgorithme passe e 
retape E50 de memorisation de I'image comgee. 
Comme il a dej^i ete dit, en variante, les etapes prin- 
cipales E10 e E40 peuvent etre derouiees plusieurs 
fols sequentiellement avant de passer d retape E50. 

En reference d la figure 8, un proc6de de 
correction de deformation d'image due d la derive 
est illustre sous la forme d'un algorithme compose de 
cinq etapes principales E100 e E500, analogues aux 
etapes E10 d E50 precedemmentdecrites en referen- 
ce d la figure 4. La principale difference est que la 
correction de rimage est effectuee colonne par co- 
lonne, et non ligne par ligne. A Uetape E100 sont lus 
N X M niveaux de luminosite de pixel memorises d'une 
image d N lignes et M cdonnes. 

L'etape E200 determine et memorise les decala- 
ges des cdonnes prtees deux d deux successrve- 
ment uetape E200 est composee d'un ensemble 
d'etapes E210 d E270 qui serontdecrites en detail en 
reference e la f^ure 9. 

L'etape E300 traite les decalages entre colonnes 



calage non lie e la prise de vue, mais existent reelle- 
ment sur la bande de terrain photographiee. L'etape 
E300 est composee d'un ensemble d'etapes E310 e 
E340 qui seront decrites en detail en reference d la fi- 
gure 10. 

L'etape E400 detenmine rimage corrig6e d partir 
des decalages entre colonnes determines e retape 
E300. En variante. l'etape E400 determine rimage 
comgee e partir des decalages entre colonnes deter- 
mines e retape E200; cette variante ne comprend pas 
retape E300. 

Dans tous les cas, retape E400 se compose d'un 
ensemble d'etapes E410 d E460 qui seront decrites 
en detail en reference e la figure 11. 

La nouvelle Image conigee des effets de la deri- 
ve est memorisee e l'etape E500etest, de preference 
egalement aff ich6e e I'ecran d'une console de visua- 
lisation associee au calculateur. 

En variante, rimage comgee d retape E400 peut 
etre traitee e nouveau e partir de retape El 00. 

Comme montree e la figure 9, la determination 
des decalages entre colonnes comprend les etapes 
E210 e E270, qui sont respectivement analogues aux 
etapes E21 d E27 decrites en reference d la figure 5, 
d la difference pres que les decalages ne sont pas 
determinees ligne e ligne, mais colonne e cdonne. 

La lettre "C" est ajoutee aux noms des differents 
parametres. Les cdonnes sont comparees deux d 
deux pour determiner un premier decalage DecCmin 
de la cdonne C„ par rapport e la colonne ^ . ex- 
prime en nombre entier de pixel, puis un second de- 
calage DC^ de la colonne par rapport d la colonne 
C,^ , expnme en nombre non entier de pbcel. Loreque 
tous les decalages DC„. pour2gm^M, sont determi- 



11 



EP 0 629 082 A1 



12 



n6s, Talgorithme passe ^ r^tape E310. 

En r6f6rence ^ la figure 10, la traitement des d6- 
calages entre coionnes comprend quatre 6tapes 
E31 0 ^ E340. Ce traitement a pour but de corriger les 
valeurs des d^calages DCm* po^^ 2^n^M, d^termi- 
n6s k r^tape E260. 

L'6tape E310 comprend le calcul de la somme SC 
et de la moyenne MC des d^calages DC^, pour 
•2^m^M. 

L'^tape E320 Initialise le param^tre m d 2 et r6ta- 
pe E330 determine le d^calage corrig6 DCCm de la co- 
lonne C^, qui est 6gal k la nrK)yenne MC des d6cala- 
ges DC„,. Le param^tre m est incr^ment^ de 1 pour 
passer d la colonne suivante. AT^tape E340, si m est 
inf^rieur ou 6gal d M, I'^tape E330 est d nouveau d6- 
roul6e pour cette valeur de m; sinon, Talgorithme pas- 
se d r^tape suivante E340. 

Cette correction de derive suppose que la derive 
de ra^ronef est quasiment constante pendant la prise 
de vue de Tlmage ce qui implique que la correction est 
identique pour toutes les coionnes. 

En variante lorsque la derive de I'a^ronef est mo- 
difiSe graduellement pendant la prise de vue de i'ima- 
ge, le calcul de la somme SC et de la moyenne MC 
des d6calages DCm* pour 2^m^M, est remplac6 par 
le calcul d'un ensemble de sommes et moyennes cal- 
cul§es sur des groupes de valeurs de dScalage. Par 
exemple, lorsque la derive pr^ente trois valeurs suc> 
cessives pendant la prise de vue, les sommes et 
moyennes sont 

ml 

SCi = X DCj, ; 
MCi = SCi/m1 

MC2 = SC2^(m2 - ml) 



= i DC^ ; 

MC3 = SC3/(M . m2) 
avec 2^m1^m2^M. Le d6calage corrig6 DCc„, est 
alors 6gal soit d SCi , soit ^ SC2. solt ^ SC3 suivant que 
m est compris entre 2 et ml , ou ml 4 1 et m2, ou m2 
+ 1 et M. 

En variante, Talgorithme ne comprend pas les 
stapes E310 d E340 pr6c6demment d^crites. Dans 
ce cas, ralgorithme passe directement de I'^tape 



E270dr6tape E410. 

En rSf^rence d la figure 11, la determination de 
I'image corrig6e comprend les stapes E410 d E460, 
respectivement analogues aux stapes E41 ^ E46 d6- 

5 crites en reference k la figure 7. La determination de 
Timage corrig6e n'est pas effectu^e ligne par ligne, 
mats colonne par colonne. Pour chaque colonne C^ 
sont determines les niveaux de luminosite corrtges de 
pixel PCn^, pour l^n^N, chacun en fonction de ni- 

10 veaux enreglstres ou calcuies ^ I'etape E42 
PIXn4ECm,m et PIX„^„^i,m de doux pixels adjacents 
dans la meme colonne. Lorsque toutes les coionnes 
sont corrlgees, Talgorithme passe d I'etape E500 de 
memorisation de Timage corrigee. 

IS Les deux types de correction d'image respective- 
ment suivant les etapes E10 d E40 et El 00 d E400 
sont par exempie effectues sequentiellement sur une 
image, suivant cetordreou I'ordre inverse, ou bien un 
pperateur choisit de n'ef fectuer que l*un des deux ty- 

20 pes de correction pour une Image donnee. 



Rdvendlcatlons 

2$ 1 - Precede de correction spatiale d'image nume- 
rique composee d'un ensemble de rangees de pixels 
(Lt d Ln; Ci d Cm) qui sont des lignes ou des coionnes 
de rimage, 

caracterise en ce qu'il comprend 

30 la comparaison (E23;E230) des rangees deux 

^ deux pour determiner (E25;E250) des premiers de- 
caiages entre rangees (Decmin;DecCmin) exprimes 
en nombre entier de pixel, 

la determination (E26;E260) de seconds deca- 

35 lages entre rangees (D„;DCm) exprimes en nombre 
non entier de pixel, en fonction des premiers decala- 
ges, 

la determination (E42;E420) pour chaque ran- 
gee (Ln;Cm) d'un troisieme decalage (Decaln;De- 

40 calC J de ladite rangee (Ln;CJ par rapport k une ran- 
gee de reference (Li;Ci) en fonction des seconds de- 
calages (DniDCJ, et 

la correction de chaque rangee (E44;E440) 
par un decalage des pixels de ladite chaque rangee 

45 par rapport d la rangee de reference egal k la partie 
entiere du troisieme decalage pour ladite chaque ran- 
gee, et par un lissage spatial des pixels de ladite cha- 
que rangee avec des coefficients de lissage 
(kni»kn2:kCmi,kCm2) fonction du troisieme decalage 

50 pour ladite chaque rangee. 

2 - Precede de correction conforme k la revendi- 
. cation 1, caracterise en ce que la comparaison 
(E23;E230) des rangees deux e deux comprend pour 
deux rangees donnees (Ln,Ln.t;Cm,C„).i), le calcul de 

55 somnries de correlation (Soeo'SDeoc) entre les deux ran- 
gees decaiees respectivement lorsque Tune est soumi- 
se par repport e I'autre e des decalages predetermines 
exprimes en nombre entier de pixel (Dec;DecC), la re- 
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cherche d'"ne sornme de correlation mini,r«le (So,. des seconds d^calages. 

mto.5>Deccmin) pamii les sommes de corr6lation calcu- 
I6es afin d'en d6duire ledit premier d6caiage (Dec- 
min;DecCmin) entre les deux rang^es exprim6 en 
nombre entier de pixel. ^ 

3 - ProcW6 de coTOCtfon conforme ^ la revendf- 
catlon 2. dans lequel une somme de correlation 
(SDec;SDocc) entre deux rang6es d6cai^ d'un d6ca- 
lage (Dec;DecC) expnm6 en nombre entierde pixel est 
6gale d la somme des valeurs absolues des differences io 
des valeurs des pixels d'une rang^e et des valeurs des 
pixels d^caies dudit d^calage (Dec;DecC) de Tautre 
rangee. 

4 - Proc6d6 de correctbn confonne d Tune quel- 
conque des revendications 2 ou 3, dans lequel la d6- is 
termination (E26;E260) du second d6calage entre 
deux rangees (D„;DCJ exprimSe en nombre non en- 

tier de pixel est r^alisee par une interpolation des 
sommes de correlation (SoecSo^^) fonction des de- 
calages (Dec;DecC) pr6detennines entre les deux 20 
rangees, par une fonction polyndmiale d'ordre supd- 
rieur ou egal A 2 autour du premier decalage (Deo- 
min;DecCmin), le second decalage (Dn;DC^) etant 
I'abscisse tfun minimum de la fonction d'interpola- 
tion. 

5 - Precede de correction conforme d Tune quel- 
conque des revendications 1 d 4. dans lequel la de- 
termination (E42;E420) pour une rangee donnee 
(Ln;Cm) d'un troisieme decalage (Decaln;DecalCJ de 
ladite rangee (L„;C^) par rapport d une rangee de re- 
ference (L,;Ci) est realises par la somme des se- 
conds decalages (D„;DCJ des rangees comprises 
entre la rangee de reference (Li;C,) et la rangee don- 
nee (U;C„). 

6 - Proc6d6 de correction confonme e Tune quel- 
conque des revendications 1 d 5. dans lequel les 
coefficients de lissage spatial iK^.K2*C^,MC^ 
d une ranges (Ln;CJ sont tels que leur somme est 
egaie d 1. 

7 - Procede de correction conforme Tune quel- 
conque des revendications 1 d 6, dans lequel les 
coefficients de lissage spatial (kn„ k^;kC„,.kC^ 
d une rangee (U,;C J sont tels qu'un premier coeffi- 
cient (k„i;kC„,i) est egal d la partie decimals du troi- 
sieme decalage (Decal„;DecalC^ de la rangee 
(U;C J et qu'un second coefficient (KaikC^^ est egal 
au compiementaire k 1 du premier coefficient. 

8 - Procede de correction conforme d Tune quel- 
conque des revendications 1 d 7, dans lequel les se- 
conds decalages entre rang6es (D„) sont remplaces 
(E33) pardes decalages corriges (DcJ qui dependent 
des seconds decalages et dont leur valeur moyenne 
est nulls. 

9 - Precede de conrection spatiale dimage 
conforme ^ Tune quelconque des revendications 1 d 55 
7, dans lequel les seconds decalages entre rangees 
(DCJ sont remplaces {E330) par des decalages 
corriges respectlvement egaux d la moyenne (MC) 
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